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We theoretically map out the ground state phase diagram of interacting dipolar fermions in one- 
dimensional lattice. Using a bosonization theory in the weak coupling limit at half filing, we show 
~1 " that one can construct a rich phase diagram by changing the angle between the lattice orientation 

and the polarization direction of the dipoles. In the strong coupling limit, at a general filing factor, 
we employ a variational approach and find that the emergence of a Wigner crystal phases. The 
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structure factor provides clear signatures of the particle ordering in the Wigner crystal phases. 



I. INTRODUCTION 



The recent experimental progress in creating degenerate cold polar atoms/molecules with large dipolar moments 
attracted considerable attention due to the rich quantum mechanical phenomena they can exhibit [l|-|8|. In contrast 
^"i. to the contact interaction, the long-range, anisotropic dipolar-dipolar interaction between dipolar molecules offers 
promising directions for exploring novel and strongly correlated many-body physics. The experimental exploration 
(DTj. of dipolar physics started with the observation of Bose-Einstein condensation of 52 Cr and 164 Dy magnetic dipolar 
atoms [H, 0|- Later, showing promising indication of creating quantum degenerate mixtures of dipolar molecules, a 
dense gas of 40 K 87 Rb and dual-species Bose-Einstein condensate of 87 Rb and 133 Cs have been realized experimentally 



fl, |f| . The realization of dipolar molecules in an optical lattice and the first creation of quantum degenerate 
^-t dipolar Fermi gas of 161 Dy have just been reported [10 ]. For molecules with permanent electric or magnetic dipole 
moments, the range of the dipole-dipole interactions can be much larger than typical optical lattice spacings. Optical 
lattices provide rich tunable ingredients such as geometry, dimensionality, and interactions so that one can engineer 



novel many-body states [ll|. These states include various superfluid states such as p x + ip y and ci-wave superfiuid 



phases, supersolid phases, vortex lattices, various Wigner crystal phases, charge-density wave and spin-density wave 



phases [l2l - l42j |. Further, cold polar molecules in optical lattices provide a platform for novel spin models and possible 



applications in quantum computing [43l |44| 



The physics of non-polar atoms with only contact interaction in optical lattices can be reasonably described by 



the Hubbard model 45] . In the Hubbard model, atom-atom interaction is approximated by an on-site interaction U. 
' However, as the dipole-dipole interaction is long-range, the experiments of polar atoms or molecules fall outside the 
^ , range of validity of the Hubbard model. A natural extension of the Hubbard model comes from including long-range, 
■ off-site interactions between the molecules. 

The one-dimensional many-body phenomena, such as the break down of Fermi liquid theory and spin-charge 
separation, can be understood in the framework of bosonization theory 4^-48|. The bosonization theory is valid 
asymptotically at small momenta and low energies. In this Letter, we study the phase diagram of two-component, 
one-dimensional lattice fermions. While we use the bosonization theory in the weak coupling limit, a variational 
approach is employed in the strong coupling limit to study the possible Wigner crystal states. Using the bosonization 
' theory at half-filling, we show that one can achieve a rich phase diagram by changing the polarization direction with 
respect to the lattice orientation. The weak coupling phase diagram includes spin-density wave, charge-density wave, 
. — i ' singlet superfluid and triplet superfluid phases. In the strong coupling limit, at smaller filling factors, we find that 
^ ! the Ion g-range interaction induces a Wigner crystal phase. The structure factor or the density-density correlation 
function which can be measured using Bragg scattering experiments provide a clear signatures of the Wigner crystal 
phase. 

The Letter is organized as follows. In section II, we discuss the effective lattice model for the dipolar fermions in 
one dimension. In section III, we present the bosonization theory for weakly interacting fermions in the presence of 
long-range off-site interaction. Assuming that the dipoles are polarized along the applied field and taking the angle 
between the lattice direction and the applied field as a free parameter, the weak coupling limit ground state phase 
diagram at half-filling is presented in section IV. The section V is devoted to discuss the effect of inter-chain coupling 
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in realistic experimental settings. In section VI, we consider the strong coupling limit and use a variational approach 
to study the possible Wigner crystal state away from half-filling. Finally in section VII, a summary is provided. 

II. THE MODEL 

We consider a system of two-component electric or magnetic dipoles confined in a one-dimensional optical lattice 
oriented along the ^-direction. The Hamiltonian operator for the fermionic atoms in optical lattice is given by 



2m dx 2 



Vv(X) + \ v ci J tfo^j(x)V>j(x)V>j.(a:)V>t(aO 

\ ( dxdx '^ai X )^l3{ X ')Vdd{x - x')lp 7 (x')lp S (x), 



(1) 



a,/3,7,<5 

where ip^ T (x)[ip rT (x)} is a fermion field operator which creates (annihilates) a Fermi atom with mass m and pseudo- 
spin a =t, I at position x. Here the pseudo-spin a refers to the two hypcrfinc states of the atom. The optical lattice 
potential provided by the counter propagating laser is Vq(x) = Vq sin 2 (fcx), with the wave amplitude Vq and wavevector 
k = 2tt/X, where A is the laser wavelength corresponding to a lattice period d = A/2. The s-wave contact interaction 
V c i — 4iTh 2 a s /m, with s-wave scattering length a s and the effective one-dimensional dipolar-dipolar interaction Vddix), 
is related to the three-dimensional dipolar-dipolar interaction, 



V dd (r) = D* l -^± (2) 

where 9d is the angle between the ID lattice in the x-direction and the dipolar moment of the atoms align along the 
applied homogeneous electric or magnetic field in the x — z plane. The strength of the dipolar-dipolar interaction is 
D 2 = d\j (47reo) and D 2 = /io^o/ (^ 7r ) f° r electric and magnetic dipoles respectively. Here eo is the electric permittivity, 
fi is the magnetic permeability, and do is the dipolar moment. For a tight one-dimensional geometry, the level spacing 
in transverse direction is much larger than the energy per particle of the axial direction x. The integration of the 
dipolar-dipolar interaction in Eq. ([2]) over the transverse direction leads to the effective one-dimensional dipolar 
interaction 17, 20] 



v u(x) = -B.112S-SM. ( 3) 

The single atomic energy eigenstates are Bloch states and localized Wannier functions are a superpositi on of B loch 
states. For a deep optical lattice with atoms trapped in the lowest vibrational states w(x) — e~ x ^ 2l ) /Vln 1 / 2 with 
/ = \fh/{muj), the field operators i\) a can be expanded as ip a = Ci a w(x — Xi). The oscillator length I is defined 
through huj = y/AEjiVo- Here w is the oscillation frequency, obtained using the harmonic approximation around the 
minima of the optical potential well at each lattice site. The recoil energy is Er = h 2 k 2 /(2m). In terms of new 



fermionic operators Ci a at lattice site i, the effective lattice Hamiltonian for the polar fermionic system reads [111 |2 



H = -t (4a c jo- + h.c) + n it n il + X! Virrii+rrii- (4) 

The parameters t — j dxw*(x — x i)[~^-§^z + Vo(x)]w(x — Xj) is the hopping matrix element between neighboring 
sites i and j, U — Airh 2 a s J dx\w(x)\ 4 /m is the on-site interaction of the two atoms at site i, and Vi r = j dxdx'\w(x — 
Xi)\ 2 Vdd(x — x')\w(x' — x r )\ 2 is the off-site interaction of two atoms at sites i and r. Apart from this "direct" like off- 
site density-density interaction term, " exchange" like spin-spin interaction term is also present for dipolar gases 49] . 
Assuming, dc electric and microwave fields in the realistic experimental setups allow one to tune the "direct" like 
interactions to be dominant (Hoj . here we neglect the spin-spin interaction term. In terms of w, Vq, and I, the 
parameters read, t = e^ Vol{2h ^huj /2, U = Airh 2 a s / (y/2irml), and V lr = -V[l + 3 cos(20 d )]/(|i - r| 3 ). Notice that 
the tunneling energy is exponentially sensitive to the laser intensity whereas the interactions are weakly sensitive. The 
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on-site interaction can be either repulsive or attractive depending on the sign of the scattering length a s . Furthermore, 
off-site interaction Vi r , can be adjusted to be positive (repulsive) or negative (attractive) by changing the direction of 
the applied field. Here i — r ^ is a discrete variable that represents the lattice points. We consider both repulsive 
and attractive regimes under the assumption that purely attractive regime is achievable in the metastable state as has 
been experimentally demonstrated for one-dimensional bosonic Cs atoms 51 1. Perhaps, the residue of small spin-spin 
interaction term restore the mechanical stability in the attractive regime. 



III. BOSONIZATION THEORY 



For asymptotic low-energy properties and for the weak coupling reg ime of the system, the continuum limit is a 
good approximation. We use the standard bosonization techniques [52] to map the Hamiltonian into the continuum 
limit by introducing continuous fermion fields Cj CT /v« — > ipLcr(x) + ?Pr<j(x) with 



^ (a,) = i e ^» e i/^W.K-rf.] i ( 5) 



The fermion operator (x) creates a Fermi atom of pseudo-spin a on the branch 77 = R, L = ±1 of the linearized 
spectrum E(k) — vp(rjk — hp), where vp — 2dt sm(kpd) is the Fermi velocity. Here R and L refer to right movers and 
left movers, respectively. The parameter a is the standard bosonization short-range distance cut-off, which is on the 
order of a lattice constant d. The Fermi wavevector is kp — im/(2d) with particle density n. In the continuum limit, 
x = jd and the length of the chain L = Nd is finite, hence we consider the limits d — > Q and the number of atoms 
TV —> 00. The discrete variable j above represents the lattice points. The fields representing particle^ = n) and 
spin(;/ = a) fluctuations are <j) v and 9 V . They satisfy the commutator [(jj^ix), 6 v (x')] — — itt /25 fi . l ,sgn(x — x'). The 
Hermitian operators U vr7 satisfy the commutator \U n<T , U v > a '] — 2S VV ' 5 acr i . Introducing the velocities v v of particle (n) 
and spin (a) sectors and Gaussian couplings K y , and following the standard procedure, the ID particle system can 
be represented by the sine-Gordon model as [53l. l56j|. 



H =y ^ 

^ 2tt 



dx 



K v (d x e v f + — (9 X ^) 2 



(6) 
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Here we use the standard notations where 



dx cos[qVs<p n (x) + Sx] cos[qV&4>o-(x)] 



v u = M(l + W2) a -(W2) a ] 1/a (7) 
K = ri^yW2 + W2l 1/2 

[l + y^/2-^/2. 
9u = 9i\\ - .92|| T.92J- 

92u = .92|| ±52-L 
9Au = .94|| ± 94 _L 

Vu = 9v/(nv F ), 

where the upper sign refers to the particle sector (n) and the lower sign refers to the spin sector (a). In standard 
bosonization language, the coupling constants and gi± with i — 1, ...4, refer to low-energy processes of the 
interaction. The coupling g\ couples two fermions on the opposite side of the Fermi surface and the particles switch 
the sides after the interactions. This process is called backward scattering or 2kp scattering. The coupling gi couples 
two fermions on the opposite sides of the Fermi surface which stay on the same side after the scattering. This process 
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is called forward scattering. Notice that the effect of g 2 is included in the first term in Eq. (6) [56]. The coupling 
between two fermions on the same side of the Fermi surface is denoted by the coupling constant g^. The subscripts 
|| and _L refer scattering between fermions with parallel spins and anti-parallel spins, respectively. The scattering 
corresponding to the coupling constants g 3 ± and g 3 \\ occurs only in the presence of the lattice. These are the well- 
known umklapp processes where the momentum is conserved up to the reciprocal lattice vector. The parameters 5 
and q control the filling factor n = N/L. In this section of the present Letter we treat the half- filling case where (5 = 
and q — 1 so that n — 1. 

In the weak coupling limit of our model in Eq. (6), all the scattering amplitudes can be presented as follows [53W55T ] . 
The amplitudes of the backward scattering are gi± = Ud+2dJ2 x V x cos(2kpx) and gi\\ = 2dJ2 x ^ cos(2kpx). Notice 
that we use V% T — > V x to represent the discrete variable \i — r\ — > dx. The amplitudes of the forward scattering are 
g2j_ = Ud — 2d^2 x V x cos{2kpx) and g 2 \\ — ~9i\\- The amplitudes of the umklapp scattering are g^± = g±± and 
93\\ = <7i||- The amplitudes of the other scattering are gn± = g 2 ±_ and g^ = g 2 \\- For the case of weak coupling, 
the velocities and the Gaussian coupling in the particle and spin sectors are v v K v — vp, v n /K n — vp — g n /7r, and 
v a /K a =v F - g a /ir. 



IV. PHASE DIAGRAM 



In the absence of the umklapp processes and in the limit gi±_ — > 0, the Hamiltonian is quadratic. In this limit, 
various correlation functions corresponding to different quantum phases can be easily calculated. These correlation 
functions show non-universal power law decay with exponents depending on Gaussian couplings K n and K a [Hfj]. 
However, in the presence of umklapp processes and the non-zero limit of one has to treat the quantum phase 
transitions by using renormalized group theoretical techniques. 

In the present section, we consider weak coupling limit at half-filling. In the weak coupling limit, the scaling 
dimension of g$ term is always higher than that of other non-linear terms in our model 53, 56|. Therefore, we set 
<?3|| = in the present study. The effect of <?ij_ and g3± is taken into account using renormalized group (RG) equations 
as has been done in Ref. [53T |. Changing the cut-off a — > e dl a with I = ItlL, the RG equations within one-loop order 
is given by 



^ = -vUD (8) 
dy V 4>(l) 



dl 



where y„o(0) = 2(K» - 1), 2/o-o(0) = 2(K a - I), y, l( /,(0) = g3j_/(yrv n ), and ^(0) = gix/{^v a ). Notice that we consider 
the weak coupling limit at half-filling. At these limits, the RG equations for particle and spin sectors are decoupled [56j . 
These equations determine the RG flow diagrams as presented in FIG. 2 of ref. [53[- The RG equations for velocities 
at these limits are trivial and velocities have no effect on scaling dimension. Following the same arguments as in ref. 
[55} . the weak coupling phase diagram at half-filling can be extracted from the RG flow diagram. 
For the spin-gap transition [y = a), Eq. (8) gives 



lyao(0) + 1 

This shows that the spin gap opens when y a o(l) < 0. For the weak coupling limit, where gu/i^vp) <C I, Gaussian 
coupling K v = [1 — g v j '(nv^)]^ 1 ^ 2 can be approximated as K v ~ I + g u /(2nv v ). In this limit, the condition y<jo{l) < 
translates into g„ < 0. Therefore, the phase boundary between the charge-density wave (CDW) and the spin-density 
wave (SDW) phases at half-filling is determined by U — 4V[1 + 3 cos(26>d)] X)m=i( — l)"V m3 - The sum over m controls 
the range of the long-range interaction. For example, m = I represents only the nearest neighbor interaction. 

On the other hand, the condition for the charge gap is g3± > \g n \. This condition gives two possible phase 
boundaries; one is g^± — —g n < and the other is <?3_l = g n > 0. As there is no continuous symmetry breaking in 
one dimension, these phase transitions, due to the opening of a charge gap, are not true phase transitions but have 
power law decaying correlations. These Bcrcnzinskii-Kostcrlitz-Thouless type transitions are due to the SU(2) and 



hidden SU(2) symmetries in the particle (charge) sector [571. l58j|. 
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FIG. 1: The phase diagram of one-dimensional polarized dipolar fermions in the weak coupling limit. The angle 0d 

is the polarized angle with respect to the ID lattice orientation. We set the filling factor to n = 1 and the 
long-ranged dipolar interaction up to 100 lattice sites. The abbreviated phases are SDW: spin-density wave, CDW: 
charge-density wave, TSF: triplet superfluid phases and SSF: singlet superfiuid. 



The boundary between CDW phase and the singlet super-conducting correlation (SSF) phase is given by the 
conditions U < and ^3 + 3cos(26>d)] cos(2kpx) = 0. At half-filling, this condition translates into U < and 
cos(20 d ) = -1/3. 

The phase boundary between SDW phase and the triplet supper-conducting correlation (TSF) phase is determined 
by the conditions U > and U — — 4V[1 + 3cos(20<2)] ^ m=1 (— l) m /m 3 . Similar to the quadratic Hamiltonian [56j j . 
a Gaussian type transition take place between two gapped phases when y n § = and y n o < 0. Since non-linear term 
vanishes on this Gaussian type transition line, this transition between SSF and TSF phase does not emerge from the 
RG equations (53|. Instead, the scaling dimensions on the Gaussian line determines the transition at g n < at g 3 ± = 0. 
Therefore, the phase boundary between SSF and TSF phases are given when U = 4V[1 + 3 cos(20d)] X)m=i(~l)"V m3 
and cos(20d) < —1/3. The resulting phase diagram in U — 0d plane for 2V — 1 is shown in FIG.Q] 

All the phases in the rich phase diagram in FIG. [T] can be constructed experimentally just by fixing the on-site 
interaction (i.e fixing the laser intensity of the counter propagating lasers) and then changing the polarization direction 
with respect to the lattice orientation. However, the interactions have to be weak and the filling factor must be unity. 
By controlling the total number of particles in experiments, the filling factor at the center of the lattice can be set 
to unity. In-situ density imaging or Bragg spectroscopy can be used to distinguish the charge-density wave and 
the spin-density wave phases. The superfluid phases can be detected via pair correlation measurements using noise 
spectroscopy [59| . 

Notice that the boundary between the singlet superfluid and charge-density wave phases does not depend on the 
range of the interaction. However, all of the other boundaries have an effect on the range of long-range dipolar 
interaction. The shift of the boundaries due to the long-range part of the interaction is shown in FIG. [2] for a fixed 
V = 0.5. The interaction strength V is always positive so that the qualitative features of the phase diagram do not 
change with V. 
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FIG. 2: The effect of the interaction range. Panel (a) shows the boundary line between a charge-density wave phase 
and a spin-density wave phase for m — 50 (black) and m — 1 (gray). Panel (b) shows the boundary line between a 
spin-density wave phase and a triplet supcrfluid phase for m — 50 (black) and m = 1 (gray). Notice that m = 1 

represents only the nearest neighbor interaction. 

V. THE EFFECT OF INTER-CHAIN COUPLING 

In the field of cold-atomic physics, one-dimensional systems are realized by creating an array of many ID-tubes. Even 
though the tunneling between tubes is absent for well separated chains, the long-range dipolar-dipolar interaction can 
still cause coupling between tubes. In the presence of inter-chain coupling, we must modify our original Hamiltonian 
in Eq. (1) by adding the inter-chain interaction term, 
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FIG. 3: The variational parameter r\ for the quarter filling case. 



Hi = - J dxdx^ljx^ ^V^x-x^^x^six), (10) 

where j ^ f is the chain index. Generalizing the ID Wannier function w(x) — > w{x,y) = e _ [ 2:2 +(2' _ J- R ) 2 ]/(2' 2 ) 
the inter-chain coupling can be approximated by V± = 2D 2 sin 9%/ R 2 , where we assume that the neighboring chain 
is R distance away in the y-direction. In the absence of the lattice, the planar array of ID tube systems has been 
studied using bosonization theory [Hoi ]. It has been shown that the inter-chain interaction is irrelevant, except when 
6d — C where the long-range interaction vanishes along the lattice in the x-direction. When 9d — 6 C , the ID system 
approaches the boundary between CDW and SSF phases and the long-range positive interaction between neighboring 
chains induces a type of CDW phase in the transverse direction 60]. Even in the presence of a lattice, the inter-chain 
coupling can induce an inter-chain CDW phase. As a result, the intra-chain CDW-SSF phase boundary shifts toward 
the CDW phase allowing the SSF phase to stabilize over the CDW phase into a larger region of the phase diagram. 



VI. THE STRONG COUPLING LIMIT 



As has been shown above, for any positive on-site and off-site interactions, the particle system produces a insulating 
phase at half-filling. For any commensurate filling factors away from half-filling, the umklapp scattering is an irrelevant 
perturbation (6lj . For the ID Hubbard model without the off-site interactions, the system remains in a metallic phase 
at any filling factor away from half-filling. However in the presence of long-range interactions, when the average 
particle spacing 1/n is comparable to the range of the interaction, atoms may form a self-organized pattern known as 
a Wigner crystal phase. The transition into this insulating phase occurs at Luttinger parameter K n = n 2 . 

For a continuous ID system, the form of l/r° interaction has been studied in bosonization theory by treating the 
long-range forward scattering as a perturbation [|32j]. This study shows that for /3 > 0, the forward scattering is an 
irrelevant perturbation. For /3 = 1, the Fourier transform of the interaction has a logarithmic divergence. At this 
limit, the bosonization theory predicts the existence of quasi- Wigner crystal phases with AUf density correlations [(32}. 
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At low filling factors, these 4kp correlations are dominant over the 2kp Friedel density correlations. The study in 
Ref. J62J is based on a perturbation approach. Therefore, the existence of quasi- Wigner crystal at larger interactions 
for /3 = 3 is not ruled out. Indeed, by investigating the structure factor using an exact diagonalization method, the 
existence of Wigner crystals at strong coupling limits and lower filling factors has been verified in Ref. [25j . This 
verification can be justified by the conditions that the dipolar-dipolar interaction Vdd{r) > 0, Vdd(r) — > as r — > oo, 
and Vdd(r + 1) + Vdd(r — 1) > 2Vdd{f) for any r > 1. Since the bosonization techniques rely on a linear band dispersion, 
low energies and long wavelengths, our approach above does not predict the Wigner crystal phase even away from 
half-filling. 

In order to study the possible existence of a Wigner crystal phase at low filling factors, we use a variational 
approach. Let's set 9d = so that the dipolar-dipolar interaction is purely repulsive. We consider the strong coupling 
regime where U, Vi r ^> t. As our motivation is to study whether Wigner crystal phases are favorable due to the 
long-range interaction, we consider the limit U — > oo. In this limit, both double occupancy and mixing of different 
spin configurations are eliminated so that we can neglect the local interaction term and suppress the spin index of 
the operators in Eq. (j4]). In other words, the Wigner crystal phase associates with charge ordering at this limit can 
be described as spinless fermions on the lattice. In Fourier space, the resulting Hamiltonian reads, 



H = ^2 e k4.c k + ^- ^2 V(q)n q n- q , (11) 

k q 

where = — 2t cos(fcd), n q ~ c], Cfc, and V(q) = 2(qd)Ki(dq) is the Fourier transform of the off-site interaction 



"k+q 

of the form Vi-j — V[L/n sin(-7r|i — j\/L)]~' A . Here K\{x) is the modified first order Bessel function [63j. The open 
boundary conditions in realistic cold atoms systems may cause edge localization phenomena if the number of lattice 
sites is small [64 |. However, in the presence of large number of lattice points, we believe that these effects are absent. 
Therefore, we assume that the optical lattice obeys periodic boundary conditions and introduced a chord distance 
between sites i and j to include the periodic boundary conditions. We consider any commensurate filling factors in 
the form n — N/L — 1/s so that one dipolar particle is occupied in a periodic sequence of unit cells of size s in the 



Wigner crystal phase. Following Ref. 65(, we take our variational wave function in the form \ip{ri)) = exp[— r/TWipo), 



where, T = —l/(t) J^ fc e^c^Cfe, 77 is the variational parameter and 



1^0)= n ^(4+4+q)|o>. (12) 

keRBZ k 

Here Q — 2w/s, |0) is the vacuum state, and RBZ stands for Reduced Brillouin Zone. This wave function is analogous 



to the well- known Gutzwiller wave function 66( which is used to explain the metal-Mott-insulator transition at half 
filling. Similar to the suppression of doubly occupied states in Gutzwiller wave function, the exponential operator 
sitting in front of our variational wave function suppresses high kinetic energy states. The normalization factor 



9 



N% = cxp[— 2-qek/t] + exp[—2rjek+Q/t} = A\ + B\. Notice that |V>o) is proportional to the classical ground state of 
the Wigner crystal phase in the absence of tunneling between sites. The variational parameter r\ is determine by 
minimizing the ground state energy E g = (ij)(r])\H\i()(r])) = (KE) + {V), where the first term is the kinetic energy and 
the second term is the off-site interaction energy. By converting the sums in to integrals over the RBZ, the variational 
kinetic and interaction energies take the form, 

The structure factor S(q) = (n q n^ q ) above has the form 



S{q) - + UJ J RBZ 2, NlNl_ q J RBZ n Nt ' (14) 

Notice that the q sum in the interaction energy run over the entire Brillouin zone (FBZ) including q = 0. This is 



different from electronic systems where the q = term is omitted due to the divergency of interaction 6jj. The 
positive background charges in the electronic lattice ensures the cancelation of this divergency. As a demonstration, 
the variational parameter rj for the quarter filling case (n = 1/2) is shown for different values of the interaction 
strengths in FIG.[3l Notice that the variational parameter reaches the classical Wigner crystal phase limit (77 = 0) for 
larger interaction strengths while it reaches the liquid phase value (77 — > 00) in the opposite limit. As our variational 
wave function always represents the Wigner crystal phase, our approach does not allow us to study the phase transition 
between Wigner crystal and liquid phases. The qualitative behavior of r\ is similar for other filling factors, however 
as the filling factor decreases, the variational parameter r\ increases. A justification of phase transition from a liquid 
phase to a Wigner crystal phase is provided at the end of this section. 

As the density-density correlation function is related to the structure factor, the evolution of the structure factor in 
FIG. (U|) shows how the particle modulation builds up as one increases the interaction. We have shown the results for 
two filling factors n = 1/2 and n — 1/4 corresponding to Q = n and Q = tt/2, respectively. The qualitative behavior 
for other filling factors are the same. The reduction of the structure factor at higher interaction is due to the transfer 
of some of its weight to the Bragg peak at q — Q. The weight transferred to the Bragg peak (I s ) can be calculated 
using I s = n — [S(Q) — n 2 ], where n = Q/(2tt). This peak intensity as a function of the interaction for the quarter 
filling case is shown in FIG. [5] As one expects, the peak intensity goes to zero for non-interacting systems. Since our 
variational approach is valid only for the Wigner crystal phase, the peak intensity is always non-zero for any finite 
interactions. However, if the Wigner crystal phase is absent, then the peak intensity must be zero. Experimentally, 
these Bragg peaks can be probed by measuring the structure factor using Bragg scattering or imaging techniques 



Since the Bragg peak at q — Q corresponding to the periodic ordering of the Wigner crystal phase, the average 
density distribution of the lattice can be written as n(x) = n + I s cos(Qx). This quantity for different interaction 
parameters is shown in FIG.|BJ As can be seen from the figure for both n = 1/2 and n = 1/4 filling factors, the higher 
interactions enhance the peak structure showing the crystalline structure in the density distribution. This periodic 
density order can be probed by using a currently available experimental technique, known as quantum gas microscopy 



As we discussed in Sec. V, for a more realistic experimental setups, one has to consider the inter-chain interaction 
in the form given in Eq. (|10[) . For 9^ — intra-tube interaction is repulsive, however fermions in different tubes 
can attract or repel depending on their dipolar moment alignment and the tube separation. For attractive inter-tube 
interactions, the system forms a clustered Wigner crystal phase j39[. This phase is coherent and Wigner crystals in 
both tubes locked to each other. On the other hand, for repulsive inter-tube interactions, the Wigner crystal phase is 
in in-coherent state. 

As we mentioned before, our variational wave function always represents a Wigner crystal phase. In order to justify 
the phase transition between a liquid state and a Wigner crystal phase, here we compare the energy of the Wigner 
crystal phase and the liquid phase perturbatively. Taking the unperturbed wave function as a free particle state, the 
liquid state energy in the first order perturbation is given by 
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FIG. 5: The weight of the Bragg peak (I s ) for the quarter filling case. 




(a) Filling factor n = 1/2 (b) Filling factor n = 1/4 



FIG. 6: [Color online] The density variation in the Wigner crystal phase for two different filling factors. 
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By comparing the liquid state energy E\ and the Wigner crystal state energy E g , we find that the liquid state is 
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favorable for small V values and phase transition between these states takes place at a finite V values for all filling 
factors. For example, for the case of quarter filling case, we find the phase transition at V = 4.74t. 

VII SUMMARY 

We have studied dipolar fermions in a one-dimensional lattice using the bosonization theory and a variational 
approach. In the weak coupling limit at half- filling, the bosonization theory predicts the appearance of several quantum 
phases as one change the polarization direction of the dipoles relative to the one-dimensional lattice orientation. The 
quantum phase diagram includes a charge-density wave, a spin-density wave, a singlet superfluid, and a triplet 
superfluid phases. In the strong coupling limit at lower filling factors, our variational method predicts the emergence 
of a Wigner crystal phase due to long-range interaction. The structure factor and the density distribution clearly 
indicates the existence of Wigner crystal at larger interactions. The entire rich phase diagram resulting from the 
competition between kinetic energy and the on-site and off-site long-range interactions can be detected by using 
currently available experimental techniques. 

VIII ACKNOWLEDGMENTS 

We thank Erik Weiler for carefully reading the manuscript. 



[1] T. Lahaye, T. Koch, B. Frohlich, M. Fattori, J. Metz, A. Griesmaier, S. Giovanazzi, and T. Pfau, Nature 448, 672 (2007). 
[2] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, D. S. 

Jin, and J. Ye, Science, 322, 231 (2008). 
[3] J. Stuhler, A. Griesmaier, T. Koch, M. Fattori, T. Pfau, S. Giovanazzi, P. Pedri, and L. Santos, Phys. Rev. Lett. 95, 

150406 (2005). 

[4] Emily Altiere, Donald P. Fahey, Michael W. Noel, Rachel J. Smith, and Thomas J. Carroll, Phys. Rev. A 84, 053431 

(2011) . 

[5] K.-K. Ni, S. Ospelkaus, D. J. Nesbitt, J. Ye, D. S. Jin, Physical Chemistry Chemical Physics 11, 9626 (2009). 
[6] D. Wang, B. Neyenhuis, M. H. G. de Miranda, K.-K. Ni, S. Ospelkaus, D. S. Jin, and J. Ye, Phys. Rev. A 81, 061404(R) 
(2010). 

[7] Mingwu Lu, Seo Ho Youn, and Benjamin L. Lev, Phys. Rev. Lett. 104, 063001 (2010). 

[8] D. J. McCarron, H. W. Cho, D. L. Jenkin, M. P. Kppinger, and S. L. Cornish, Phys. Rev. A 84, 011603(R) (2011). 
[9] Amodsen Chotia, Brian Neyenhuis, Steven A. Moses, Bo Yan, Jacob P. Covey, Michael Foss-Feig, Ana Maria Rey, Deborah 
S. Jin, and JunYe, Phys. Rev. Lett. 108, 080405 (2012). 
[10] Mingwu Lu, Nathaniel Q. Burdick, and Benjamin L. Lev, PRL 108, 215301 (2012). 

[11] For tutorials on ultracold dipolar gases in optical lattices, see for example; C Trefzger, C Menotti, B Capogrosso-Sansone 
and M Lewenstein, J. Phys. B: At. Mol. Opt. Phys. 44 (2011) 193001 AND M.A. Baranov, Physics Reports 464 (2008) 
71. 

[12] K. Goral, L. Santos, and M. Lewenstein, Phys. Rev. Lett. 88, 170406 (2002). 
[13] G. Pupillo et al., Phys. Rev. Lett. 100, 050402 (2008). 

[14] D.W. Wang, M.D. Lukin, and E. Dernier, Phys. Rev. Lett. 97, 180413 (2006). 

[15] C. Trefzger, C. Menotti, and M. Lewenstein, Phys. Rev. Lett. 103, 035304 (2009). 

[16] S. Zollner, G. M. Bruun, C. J. Pethick, and S. M. Reimann, Phys. Rev. Lett. 107, 035301 (2011). 

[17] F. Deuretzbacher, J. C. Cremon, and S. M. Reimann, Phys. Rev. A 81, 063616 (2010). 

[18] R. Citro, E. Orignac, S. DePalo, and M. L. Chiofalo, Phys. Rev. A 75, 051602(R) (2007). 

[19] A. S. Arkhipov, G. E. Astrakharchik, A. V. Belikov, and Y. E. Lozovik, JETP Lett. 82, 39 (2005). 

[20] S. Sinha and L. Santos, Phys. Rev. Lett. 99, 140406 (2007). 

[21] R.-Z. Qiu, S.-P. Kou, Z.-X. Hu, X. Wan, and S. Yi, Phys. Rev. A 83, 063633 (2011). 
[22] S. Yi, T. Li, and C. P. Sun, Phys. Rev. Lett. 98, 260405 (2007). 

[23] A. E. Golomedov, G. E. Astrakharchik, and Yu. E. Lozovik, Phys. Rev. A 84, 033615 (2011). 

[24] B. Capogrosso-Sansone, C. Trefzger, M. Lewenstein, P. Zoller, and G. Pupillo, Phys. Rev. Lett. 104, 125301 (2010). 
[25] Zhihao Xu and Shu Chen, Phys. Rev. A 85, 033606 (2012). 

[26] S. G. Bhongale, L. Mathey, Shan- Wen Tsai, Charles W. Clark, Erhai Zhao, Phys. Rev. Lett. 108, 145301 (2012). 

[27] Tomasz Sowinski, Omjyoti Dutta, Philipp Hauke, Luca Tagliacozzo, and Maciej Lewenstein, Phys. Rev. Lett. 108, 115301 

(2012) . 

[28] Liang He and Walter Hofstetter, Phys. Rev. A 83, 053629 (2011). 
[29] K. Mikelsons and J. K. Freericks, Phys. Rev. A, 83, 043609 (2011). 



12 



[30] Chungwei Lin, Erhai Zhao, and W. Vincent Liu, Phys.Rev.B. 81. 045115 (2010); Phys.Rev.B. 83 119901 (2011). 

[31] M. M. Parish and F. M. Marchetti, Phys. Rev. Lett. 108, 145304 (2012). 

[32] F. M. Marchetti, M. M. Parish, pre-print laxXiv: 1207.40681 

[33] N. R. Cooper and G. V. Shlyapnikov, Phys. Rev. Lett. 103, 155302 (2009). 

[34] A. Pikovski, M. Klawunn, G. V. Shlyapnikov, and L. Santos, Phys. Rev. Lett. 105, 215302 (2010). 
[35] N. T. Zinner, B. Wunsch, D. Pekker, and D.-W. Wang, Phys. Rev. A 85, 013603 (2012). 
[36] M. Klawunn, J. Duhme, and L. Santos, Phys. Rev. A 81, 013604 (2010). 

[37] B. Wunsch, N. T. Zinner, I. B. Mekhov, S.-J. Huang, D.-W. Wang, and E. Dernier, Phys. Rev. Lett. 107, 073201 (2011). 
[38] P. Lecheminant and H. Nonne, Phys. Rev. B 85, 195121 (2012). 

[39] Michael Knap, Erez Berg, Martin Ganahl, and Eugene Dernier, Phys. Rev. B 86, 064501 (2012). 

[40] Kai Sun, Congjun Wu, and S. Das Sarma, Phys. Rev. B 82, 075105 (2010). 

[41] Y. Yamaguchi, T. Sogo, T. Ito, T. Miyakawa, Phys. Rev. A 82, 013643 (2010). 

[42] N. T. Zinner, G. M. Bruun, Eur. Phys. J. D 65, 133 (2011). 

[43] A. Micheli, G.K. Brennen, P. Zoller, Nature Physics, 2, 341-347 (2006). 

[44] D. DeMille, Phys. Rev. Lett. 88, 067901 (2002). 

[45] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[46] Quantum Physics in One Dimension, T. Giamarchi, Oxford University Press, 2004. 

[47] V. J. Emery, in Highly Conducting One Dimensional Solids, edited by J. T. Devreese et al. (Plenum, New York, 1979), p. 
327. 

[48] M. A. Cazalilla, J. Phys. B: At. Mol. Opt. Phys. 37, SI (2004). 

[49] Kaden R. A. Hazzard, Alexey V. Gorshkov, and Ana Maria Rey, Phys. Rev. A 84, 033608 (2011). 

[50] Alexey V. Gorshkov, Salvatore R. Manmana, Gang Chen, Jun Ye, Eugene Demler, Mikhail D. Lukin, and Ana Maria Rey, 

Phys. Rev. Lett. 107, 115301 (2011). 
[51] Elmar Haller, Mattias Gustavsson, Manfred J. Mark, Johann G. Danzl, Russell Hart, Guido Pupillo, Hanns-Christoph 

Nagerl, Science 325, 1224 (2009). 
[52] F. D. M. Haldane, J. Phys. C 14, 2585 (1981). 
[53] Masaaki Nakamura, Phys. Rev. B 61, 16377 (2000). 
[54] M. Tsuchiizu and A. Furusaki, Phys. Rev. Lett. 88, 056402 (2002). 
[55] S. Capponi, D. Poilblanc, and T. Giamarchi, Phys. Rev. B 61, 13410 (2000). 
[56] Johannes Voit, Phys. Rev. B 45, 4027 (1992). 
[57] C. N. Yang and S. C. Zhang, Mod. Phys. Lett. B 4, 759 (1990). 
[58] T. Giamarchi and H. J. Schulz, Phys. Rev. B 39, 4620 (1989). 
[59] E. Altman, E. Demler, and M. D. Lukin, Phys. Rev. A 70, 013603 (2004). 
[60] Yi-Ping Huang and Daw- Wei Wang, Phys. Rev. A 80, 053610 (2009). 
[61] M. Brech, J. Voit and H. Buttner, Europhys. Lett., 12, 289 (1990). 
[62] Yasumasa Tsukamoto and Norio Kawakami, J. Phys. Soc. Jpn. 69, 149 (2000). 
[63] Hitoshi Inoue and Kiyohide Nomura, J. Phys. A: Math. Gen. 39, 2161 (2006). 
[64] Ricardo A. Pinto, Masudul Haque, Sergej Flach, Phys. Rev. A 79, 052118 (2009). 
[65] B. Valenzuela, S. Fratini, and D. Baeriswyl, Phys. Rev. B 68, 045112 (2003). 
[66] M. C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963). 
[67] H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

[68] M. Weidemuller, A. Hemmerich, A. Gorlitz, T. Esslinger and T. W. Hansch, Phys. Rev. Lett. 75, 4583 (1995). 

[69] G. Raithel, G. Birkl , A. Kastberg, W. D. Phillips and S. L. Rolston, Phys. Rev. Lett. 78, 630 (1997). 

[70] H. Miyake, et al. Phys. Rev. Lett. 107, 175302 (2011). 

[71] F. Gerbier, et al. Phys. Rev. Lett. 95, 050404 (2005). 

[72] S. Foiling, et al. Nature 434, 481 (2005). 

[73] C-L. Hung et al. New. J. Phys. 13, 075019 (2011). 

[74] Nathan Gemelke, Xibo Zhang, Chen-Lung Hung and Cheng Chin, Nature 460, 995 (2009). 
[75] Waseem S. Bakr, Jonathon I. Gillen, Amy Peng, Simon Foiling and Markus Greiner, Nature 462, 74 (2009). 
[76] Jacob F. Sherson, Christof Weitenberg, Manuel Endres, Marc Cheneau, Immanuel Bloch and Stefan Kuhr, Nature 467, 
68 (2010). 



